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1 Introduction 



The dynamical description of strongly interacting systems out of equilibrium nowadays 
is dominantly based on transport theories and efficient numerical recipies have been set 
up for the solution of the coupled channel transport equations [|I], EJ (and Refs. therein). 
These transport approaches have been derived either from the Kadanoff-Baym equations 
|§ in Refs. j|, |], H H or fr° m the hierarchy of connected equal-time Green functions 
H, |K| in Refs. |]TT] , |12|j by applying a Wigner transformation and restricting to first order 
in the derivatives of the phase-space variables (X, P) . Whereas theoretical formulations 
of off-shell quantum transport have been limited to the formal level for a couple of years 
|| |7|, [13], [14| only recently a tractable semiclassic form has been derived for testparticles 



in the eight dimensional phase space of a particle [15|, IB], [17 . 

Whereas in Refs. [15|, |IJ we have investigated the off-shell transport approach with 
respect to nucleus-nucleus collisions at GANIL, SIS and AGS energies, we here concen- 
trate on equilibration phenomena relative to the on-shell dynamics by imposing periodic 
boundary conditions for the system confined to a box of size V = L 3 , where L denotes 
the length of the cubic box. We, furthermore, compare the equilibrium nucleon and A 
distribution functions (t — > oo) to the statistical model (SM) employing the same spectral 
functions as in the transport approach. For related studies at higher bombarding ener- 
gies or energy densities within on-shell transport approaches we refer the reader to Refs. 

ra, pi pa, etl m 



2 Extended semiclassical transport equations 



We briefly recall the basic equations for Green functions and particle self energies that 
are exploited in the derivation of off-shell transport equations in the semiclassical limit. 

The general starting point for the derivation of a transport equation for particles with 
a finite and dynamical width are the Dyson- Schwinger equations for the retarded and 
advanced Green functions S ret , S adv and for the non-ordered Green functions S K and S > 
||. In the case of scalar bosons - which is considered in the following for simplicity - 
these Green functions are defined by 
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They depend on the space-time coordinates x, y as indicated by the indices - xy . 
functions are determined via Dyson-Schwinger equations by the retarded and advanced 
self energies E ret , Y, adv and the collisional self energy S < : 
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where Eq. (|]) is the well-known Kadanoff-Baym equation. Here denotes the (nega- 
tive) Klein-Gordon differential operator which is given for bosonic field quanta of (bare) 
mass M by = —(d^d^ + Mq); 5 xy represents the four-dimensional (^-distribution 



5 xy = 5^(x — y) and the symbol indicates an integration (from — oo to oo) as well as 
a summation over all discrete intermediate variables (cf. [?J |15|). 



2.1 The semiclassical limit 

For the derivation of a semiclassical transport equation one now changes from a pure 
space-time formulation into the Wigner-representation. The theory is then formulated 
in terms of the center-of-mass variable X = (x + y)/2 and the momentum P, which is 
introduced by Fourier-transformation with respect to the relative space-time coordinate 
(x — y). In any semiclassical transport theory one, furthermore, keeps only contributions 
up to the first order in the space-time gradients. After carrying-out these two steps the 
Dyson-Schwinger equations (@)-((D become 
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where the operator O is defined as J7|, |T3J 
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It is a four-dimensional generalization of the well-known Poisson-bracket. Starting from 
(|5|) and (^) one obtains algebraic relations between the real and the imaginary part of the 
retarded Green functions. On the other hand eq. ([F]) leads to a 'transport equation' for 
the Green function S K [15| . 

To this aim one separates all retarded and advanced quantities - Green functions and 
self energies - into real and imaginary parts, 
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The imaginary part of the retarded propagator is given (up to a factor 2) by the normalized 
spectral function 
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while the imaginary part of the self energy corresponds to half the particle width Txp- By 
separating the complex equations (|5|) and (|6]) into their real and imaginary contributions 
we obtain an algebraic equation between the real and the imaginary part of S ret , 

P 2 _ y\/f 2 _ PpVret 

RcSTp = ^ A XP . (11) 

1 XP 

In addition we gain an algebraic solution for the spectral function as 

XP ' (P 2 - M 2 - Re^pf + T 2 XP /A ' 1 1 

while the real part of the retarded propagator is given by 
p 2 _ M 2 — RpY ret 

Keb XP ~ (p2 _ M 2 _ ReT ret p) 2 + r ^ /4 " l^J 

Furthermore, the (Wigner-transformed) Kadanoff-Baym equation (|?p allows for the con- 
struction of a transport equation for the Green function S < . When separating the real and 
the imaginary contribution of this equation we find i) a generalized transport equation, 

0{P 2 - Ml - JfeES* } { S< P } - O { £< P } {ReS r x e P } 

= 2 [ ^xp ^xp ~ ^ip &XP ] > (14) 



and ii) a generalized mass-shell constraint 
[P 2 - Ml - ReE&] S< P - E< P ReST P 

= \<>{Z XP }{A XP } - ±0{T XP }{S< P }. (15) 

Besides the drift term (i.e. 0{P 2 - M 2 }{5 < } = -P^d^S<) and the Vlasov term (i.e. 
-0{PeS re *}{5<}) a third contribution appears on the l.h.s. of© (i.e. -0{£<}{PeS re *}), 
which vanishes in the quasiparticle limit and incorporates - as shown in |L5], [16], [L7| - the 
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off-shell behaviour in the particle propagation which has been neglected so far in trans- 
port studies^. The r.h.s. of flH|) consists of a collision term with its characteristic gain 
(~ T I < S > ) and loss (~ S > 5' < ) structure, where scattering and decay processes of particles 
are described. 

Within the specific term (0{T l < }{ReS ret }) a further modification is necessary. Ac- 
cording to Botermans and Malfliet |J the collisional self energy E < should be replaced by 
S < ■ Y/A to gain a consistent first order gradient expansion scheme. The replacement is 
allowed since the difference between these two expressions (S < — S < • Y/A) can be shown 
to be of first order in the space-time gradients itself Furthermore, this substitution 



is required to get rid of the inequivalence between the general transport equation (|T| 
and the general mass shell constraint flT5|) pSfl. 



Finally, the general transport equation (in first order gradient expansion) reads (15 
0,0 
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Its formal structure is fixed by the approximations applied, however, its physical contents 
is fully determined by the different self energies, i.e. -ReS^p, Yxp, ^xp that have to be 
specified in addition. 

In order to obtain an approximate solution to the transport equation (|T6|) a testparticle 
ansatz is used for the Green function S < , more specifically for the real and positive 
semidefinite quantity Fxp = AxpNxp = i S X p, 

N 

Fxp ~ £ 5®(X - X^t)) S®(P - 3(f)) S(P - e*(f)) ■ (17) 
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In the most general case (where the self energies depend on four-momentum P, time 
t = Xq and the spatial coordinates X) the equations of motion for the testparticles read 
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This also holds true for the recent numerical off-shell simulations in Refs. 
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where the notation Fm implies that the function is taken at the coordinates of the test- 
particle, i.e. F {i) = F(t, Xi(t), Pi(t), €i(t)). 

In ( P^pDD a common multiplication factor (1 — C^)^ 1 appears, which contains the 
energy derivatives of the retarded self energy 
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which yields a shift of the system time t to the 'eigentime' of particle i defined by 
ii = t/(l — C(j)). The derivatives with respect to the 'eigentime', i.e. dXi/dti, dPi/dtt 
and dei/dti then emerge without this renormalization factor for each testparticle i when 
neglecting the explicit time dependence of in line with the semiclassical approxima- 
tion scheme. The role and the importance of this correction factors have been studied in 



Ref. [16 1 for a four- momentum-dependent 'trial' potential and we refer the reader to the 
latter analysis for more details. 

Following Ref. |15|] we take M 2 = P 2 — ReY7 et as an independent variable instead of 
the energy P . Eq. (p0|) then turns to 



dM[ 
~df 



M 2 - M 2 dT (l) 



dt 



(22) 



for the time evolution of the test-particle i in the invariant mass squared as derived in 
Refs. |15|, [16| . We mention that corresponding equations of motion have been derived by 
Leupold in the nonrelativistic limit [T7J . 



2.2 Generalized collision terms for bosons and fermions 

The collision term of the Kadanoff-Baym equation can only be worked out in more detail 
by giving explicit approximations for and S > . We recall the formulation and result 
from Ref. [16|] that is based on Dirac-Brueckner theory, i.e. 



I coll (X, P, M 2 ) = Tr 2 Tr 3 Tr 4 A{X, P, M 2 )A(X, P 2 , M 2 )A(X, P 3 , M 2 )A(X, P 4 , M 4 2 ) 
\T{(P, M 2 ) + (P 2 , M 2 2 ) -> (P 3) Ml) + (P 4 , Ml))\\ s 6^(P + P 2 - P 3 - P 4 ) 
yN x p 3M 2 N x p 4M 2 f X pM 2 fxP 2 M$ ~ N x p M 2 N x p 2M 2 fxP 3 Ml /xp 4 m| ] (23) 

with 

fxPAP = 1 + vN x p M2 (24) 

and rj = ±1 for bosons/fermions, respectively. The indices A,S stand for the antisym- 
metric/symmetric matrix element of the in-medium off-shell scattering amplitude T in 
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case of fermions/bosons. In Eq. (|23| ) the trace over particles 2,3,4 reads explicitly for 
fermions 

r r ,^ ' /ft^=, (25) 

^ 2 ,r 2 (2vr) 4 J 2 y/Pi + M 2 
where a 2 , r 2 denote the spin and isospin of particle 2. In case of bosons we have instead 

Tr2 = V -±- A \ d 3 P 2 ^ } (26) 

since here the spectral function As is normalized as 
r HP 2 

J^A B (X,P) = 1 (27) 

whereas for fermions we have 
r HP 

J^A F (X,P) = 1. (28) 

We mention that the spectral function Ap in case of fermions in (|23"1) is obtained by 
considering only particles of positive energy and assuming the spectral function to be 
identical for spin 'up' and 'down' states (cf. Ref. |16|). 

Neglecting the 'gain-term' in Eq. (|2"3"| ) one recognizes that the collisional width T co u 
of the particle in the rest frame is given by 

T coll (X, P, M 2 ) = Tr 2 Tr 3 T n |T((P, M 2 ) + (P 2 , M 2 2 ) - (P 3 , M 3 2 ) + (P 4 , M 4 2 ))|\ 5 (29) 
A(X, P 2 , M 2 )A(X, P 3 , M 2 3 )A(X, P 4 , Ml) 5\P + P 2 - P 3 - P 4 ) iV^^ / X ^ M| , 

where as in Eq. (|23|) local off-shell transition amplitudes enter for the transitions P+P 2 — > 
P 3 + P 4 . We note that the extension of Eq. (|23|) to inelastic scattering processes (e.g. 
iViV — > iVA) or (7rA^ — > A etc.) is straightforward when exchanging the elastic transition 
amplitude T by the corresponding inelastic one and taking care of Pauli-blo eking or Bose- 
enhancement for the particles in the final state. The relation of the quantity Txp to the 
collisional width T co a is given by Txp = 2Pq(T deca y (X P) + T co u(XP)), where the particle 
decay width r 'decay i n the medium might also be different compared to the vacuum. 

Thus the transport approach determines the particle spectral function dynamically 
via (|29|) - with respect to the collisional width T co u - for all hadrons if the in-medium 
transition amplitudes T are known in their full off-shell dependence. Since this information 
is not available for configurations of hot and dense matter, a couple of assumptions and 
numerical approximation schemes have to be invoked in actual applications so far. 



As in Refs. Rig , [Tq| the following dynamical calculations are based on the conventional 
HSD transport approach |fj, ^] - in which PeS^p is specified for the hadrons - however, 
the equations of motion for the testparticles are extended to (p!8| , pj|j22]) . For further details 



on the elastic and inelastic transition rates we refer the reader to Ref. UTS . 
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2.3 Comment on particle number conservation 

In previous derivations of the off-shell transport equations one has started from a formula- 
tion of the non-equilibrium theory in space-time representation (x,x') and then changed 
into a phase-space representation via Fourier transformation with respect to (x — x') 
0, [T5|, [16], [I7| . The semiclassical limit then has been introduced by assuming gradients 
in X = (x + x')/2 to be small [|17[] for ReY7 xp and S%p- Here we argue that for reasons 
of symmetry in phase-space also the four-momentum derivatives in P have to be small 
to achieve a proper semiclassical limit as inherent in the formulation of the transport 
equation ( |T6"D in terms of the generalized Poisson-bracket (^). 

We briefly demonstrate in the following lines that instead of a coordinate-space repre- 
sentation one may formulate the theory equivalently in momentum-space and then trans- 
form to phase space: The momentum-dependent (two-point) functions are given as 

F PiPil = J^jdkv e^A-*V^) F w (30) 

The evolution of the retarded and advanced Green functions S ret , S adv and the Green 
function S < turns to 

(27T) 

(32) 
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with the 'kinetic' operator in momentum-space = (pf — Mq) in the case of relativistic 
scalar bosons with (bare) mass Mo. Obviously, the equations in momentum-space are 
formally equivalent to the equations in coordinate-space, i.e. the convolution integrals 
in the coordinate X2 are replaced by convolution integrals in the momentum p 2 . When 
transforming from momentum- to phase-space via a Fourier-transformation with respect 
to the four- momentum coordinate (pi — py ) 

~ J — (2tt) 4 PlPl ' ^ 

(with P = {jpi + py)/2) we gain again the familiar equations in Wigner-representation: 

q-1 nret,adv -i . — iO spret,adv nret,adv (1A\ 
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with the 'kinetic' operator in phase-space S 0XP = (P 2 + iP^d* — \d*d x — M5 2 ). Here 
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we have used the following relation for the convolution integrals in momentum-space 
f d\ Pl -p v ) iX ( Pl - Pl ,) f d 4 p 2 l0 ( s 

J (2tt) 4 7 (2^)4 *l'PiP2 **>P2Pi> ~ e *1,PX *2,PX, W 

with the Poisson-operator O defined in (|j). 

The semiclassical limit - along the conventional line of arguments - now is achieved 
by assuming gradients of the self energies Yl ret,adv in P to be small. This assumption is 
especially well taken for systems with dominantly short range interactions since the dif- 
ferent self energies T, ret,adv (p,p') become smooth functions in momentum-space such that 
a restriction to first order gradients in the momentum P = (p + p')/2 can be more easily 
justified. Note again that the four-dimensional Poisson-bracket (El) is symmetric in the 
space-time and the four-momentum derivatives. Thus a phase-space gradient expansion 
requires all gradients to be small contrary to the assumptions made in Ref. f[7 |. 



As a consequence mixed gradients of second order as d 2 /(dtdpo) have to be neglected 
in a consistent truncation scheme of first order. As shown by Botermans and Malfliet in 
the review || the particle number conservation then holds strictly. Only when keeping 
special terms of second order (~ d 2 / (dtdpo) unphysical pecularities may appear, that lead 
to a violation of particle number conservation in contradiction to the full theory (2)-(4). 



3 Infinite nuclear matter problems 

In case of infinite nuclear matter problems the solution of the transport equations simpli- 
fies since all spatial gradients (with respect to X) vanish. This implies that the momentum 
coordinates of the testparticles Pi are constant according to (19) in between collisions. 

The initial conditions of the problem then are fully specified by an initial bombarding 
energy per nucleon, that defines the relative shift in momentum of two Fermi spheres with 
a Fermi momentum Pp = 260 MeV/c, the number of nucleons N and the total volume 
of the box V = L 3 . Alternatively, one can also characterize the system by an average 



density p = N/V and energy density e (cf. Ref. Hl8| ). We use L = 10 fm and assume an 
equal number of neutrons and protons. 

In general (for t — > oo) the stationary solution of Eq. ( p3|) for a fermion h is given by 

F h (X,P,M>) = f/f'^ (37) 



with E = y P 2 + M 2 and rj = ±1 for bosons/fermions, respectively, while denotes the 
spectral function (|12"D, T the temperature of the system and \ih the chemical potential 



for the hadron. This situation corresponds to the grand canonical ensemble of quantum 
statistical mechanics. In case of infinite nuclear matter problems the dependence on X 
vanishes additionally. 
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3.1 Numerical results 



As an example for equilibration phenomena we show in Fig. [I] the time evolution of the 
quadrupole moment (involving all baryons B) 

<&>(*) = £ 7S1 / d 3 X f d*P ^ (2P 2 2 - Pi - P;) F B (X, t, P, M 2 ), (38) 

b y Z7T ) J J 2yP 2 + M 2 

for p = po and initial bombarding energies of 0.1 A GeV and 1 A GeV, respectively. In 
(|38"D the degeneracy factor is gs = 4 for nucleons and Qb— 16 for A-resonances. The solid 
lines result from the off-shell calculation while the dashed lines result from the on-shell 
limit. As can be seen from Fig. [I] the off-shell results are practically identical to the 
on-shell limits except for the very long time behaviour, where the off-shell limit needs 
some more time to achieve equilibrium. However, when defining an equilibration time r 
by a drop of the observable by the factor e -1 we find no sizeable effect from the off-shell 
propagation on r. 



The equilibrium distributions in the nucleon transverse mass M t = yPt+ M 2 are 
shown in Fig. Q for initial bombarding energies of 0.1 A GeV (upper part) and 1 A 
GeV (lower part). Again the solid lines denote the off-shell results from the transport 
calculations while the on-shell spectra are displayed in terms of dashed lines. Both limits 
(within the statistics) give the same temperature T « 97 MeV for 1.0 A GeV and T « 26 
MeV for 0.1 A GeV as can be extracted from the high energy tail of the transverse mass 
spectra. Differences can only be found for < Mo since the finite width in the nucleon 
spectral function only can show up in the off-shell case. This component is quite small at 
0.1 A GeV since the collisional width of nucleons at density po and temperature T ps 26 
MeV is about 8 MeV. 

Without explicit display we mention that the time evolution of the nucleon spectral 
function in the invariant mass M becomes broad in the initial nonequilibrium phase of the 
reaction and approaches the equilibrium distribution ( j37| ) roughly within the equilibration 
times from Fig. [1]. Since the width of the nucleon spectral function in equilibrium at 0.1 
A GeV (GANIL energy) is rather small (r « 8 MeV) we concentrate on the energy of 
1 A GeV (SIS energy) in the following, where the nucleon collisional width T co u at a 
temperature of 97 MeV is about 40 MeV and roughly 20% of the baryons are excited A 
resonances^]. 

As demonstrated above, for the equilibrated system we can extract a temperature T 
by fitting the particle spectra with the Bolzmann distribution 

^ ~ eM-Ei/n (39) 



2 Note, that in collisions of finite systems at this bombarding energy the A-abundancy is slightly lower 
due to a rapid expansion of the system and the additional compressional energy stored in the system in 
the high density phase. 
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where Ei = yp 2 + mf is the energy of particle i. We note that at the temperatures of 
interest here, the Bose and Fermi distributions are practically identical to a Boltzmann 
distribution. We find that in equilibrium the spectra of all particles (nucleons, A's and 
pions) can be characterized by a single temperature T (cf. e.g. Ref. JT^]). 



3.2 Comparison to the statistical model 

In order to investigate the equilibrium behaviour of hadron matter we compare our trans- 
port (box) calculations with a simple Statistical Model (SM) for an Ideal Hadron Gas 
where the system is described by a grand canonical ensemble of non-interacting fermions 
and bosons in equilibrium at temperature T. All baryon and meson species considered in 
the transport model (N, A, 7i~) also are included in the statistical model. 

We recall that in the SM particle multiplicities rii and energy densities e$ for particles 
with spectral functions A4 are given by 



CO 



9i [dM 7 Aj(M)/Aixp 2 dp 
Hi (27rhf J 2tt . (exp [(Ei-BiUB-SifMs) ffl-ri' 







dM J Ai(M) 4irEip 2 dp 



(2tt^) 3 J 2tt J exp [(E { - B if i B - S if i s )/T] - rj ' 







where Ei = \/p 2 + M 2 is the energy of particle i, Bi is the baryon charge, Si is the 
strangeness, and Qi is the spin-isospin degeneracy factor, while rj = ±1 for bosons/fermions, 
respectively. In Eqs. (|40|) , (fill) A*s an d A*5 are the baryon and strangeness chemical po- 
tentials. The energy density e, baryon density p and strange particle density of the whole 
system in equilibrium then is given as 

£ = $>(7> B ,//s) (42) 

i 

p = ^Bi Ui(T, y. B , ps) (43) 

i 

Ps = J2 S i n i( T ^B,Ps) = 0. (44) 

i 

As 'input' for the SM we use the same e,p and ps as in the box calculations and we 
obtain the thermodynamical parameters - T, p B , ps - by solving the system of nonlinear 
equations (P),(PD and (pi). 

We now turn to a comparison of the equilibrium distributions in mass for nucleons 
and A's, i.e. dN^/dM and dN&/dM, from the box calculations with those from the SM, 
which are obtained by integration (summation) over momentum. We first discuss the 
on-shell limit where the nucleon spectral function is represented by a 5-function around 
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the bare mass while the A spectral function - as implemented in the transport approach 
- is given by (in the A rest frame) 

A (M\ 2M 2 r 7rJV (M) 
M ] (M 2 - Ml )* + M*Tl N (M) { ' 



with 

r 7TN (M) = i,(^rU^— T f. (16) 



"q R q 2 + S 2 



where q R denotes the pion momentum in the A rest frame of the resonance and q is the 
corresponding pion three-momentum at invariant mass M, i.e. 



2 (M 2 — (M N + M 7r ) i )(M 2 - (M N - M v ) 



2 = ^ , ^ /_v_I ^ l^Z. (47) 

y 4M 2 1 J 

The quantity 5 in fl46|) is fixed by |p7 



r 2 

<5 2 = (M — Mjy — M^) 2 + (48) 

with T = 120 MeV to achieve a good description for the ttN — > A reaction in vacuum. 
Recall that in this case we have T X p/C2M) = T decay = T^jy in the A rest frame, i.e. T co u 
= 0. 

In the upper part of Fig. [3] we compare the asymptotic (t — > oo) distributions for 
nucleons (N) and A's in the on-shell limit from the transport (box) calculation (solid 
histograms) with the result from the SM at a temperature T = 97 MeV employing the 
A spectral function (f45|). Since the differences are hardly visible, we conclude that the 
transport calculation reproduces the result from the SM, where the thermodynamical 
parameters are determined by energy and baryon number conservation. We mention that 
we have discarded strangeness in this comparison, since kaons and hyperons are very 
scarce at this energy and have been switched off in both models. Since the A width ( ^6|) 
is zero below the 7rN threshold, the A mass distribution only can extend above M n + M^. 

In case of the off-shell calculation we obtain somewhat different mass distributions 
for nucleons and A's from the box calculation, which are given in terms of the solid 
histograms in the lower part of Fig. Here the nucleon spectral function becomes 
very broad and also the A distribution extends below the ttN threshold in vacuum. The 
nucleon spectral function is broadened due to the elastic (NN — > NN) and inelastic 
(NN — > NA) scattering processes, which can roughly be described by a collisional width 
T co u of 40 MeV in the nucleon spectral function 



2M 2 r 



A »w - (M2 _^Tm^ l , - (49) 
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as shown by the left dashed line in the lower part of Fig. [3]. The collisional width in case 
of nucleons is related to Txp as 2P Q T co u = Txp which reduces in the rest system of the 
particle to 2MT co u = Txp- 

The A mass spectrum in this case is more difficult to interprete. This is due to the 
fact that in the decay A — > ttN* the A may decay to a pion and an off-shell nucleon N* 
according to the nucleon equilibrium spectral function An(M), which essentially shifts 
the 7rN* threshold accordingly. However, taking into account this change in width due to 
the in-medium 7riV* decay, only, the low mass part of the A distribution is underestimated 
considerably. Here the 'collisional' channels AiV — > AN and AiV — > NN are much more 
important. We mention that the latter reaction is described by the extended detailed 



balance relation of Ref. |28|] while the elastic differential cross section is taken the same 
as for NN scattering (as a function of the momentum transfer). Since especially the AiV 
absorption reaction depends on mass M and momentum p of the baryons explicitly, it is 
not straightforward to present analytical formulas since final state Pauli blocking for the 
nucleons leads to a highly nonlinear problem. We thus have extracted the A collisional 
width r^nip, M) from the transport calculation explicitly by calculating the number of 
AiV unblocked collisions per unit time as a function of the A momentum p and mass M. 
Using (p — \P\) 



J- XP 
~2~K 



r tot ( P ,M) = r wN *(M) + rt n (p,M), (50) 



where T n N* denotes the A in-medium 7riV* width averaged over the nucleon equilibrium 
distribution function (cf. lower part of Fig. 3) and integrating over momentum with the 
appropriate Fermi function (for T = 97 MeV) we obtain the right dashed line in the lower 
part of Fig. |3| that describes the box A- distribution rather well. Thus the transport cal- 
culation at finite density po an d temperature T is consistent with the SM when employing 
the same physical baryon spectral functions. One might have expected this equivalence in 
equilibrium due to energy and baryon number conservation, however, the actual numer- 
ical result then may be regarded as a consistency test of the numerical implementation 
schemes for the various elastic and inelastic reaction channel in the medium for off-shell 
particle propagation. 

The mass integrated number of A's in the off-shell case is larger by a few % as compared 
to the on-shell case in equilibrium due to the low mass tail of the A distribution, however, 
the high mass tails of the A's are identical within the statistics achieved as demonstrated 
in Fig. [|. On the other hand, the number of pions is practically identical in both cases 
since the low mass tail of the A's is dominantly due to AiV reactions and only to a lower 
extent to the 7riV* decay as discussed above. 
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4 Summary 



In this work we have employed the semiclassical off-shell approach from Refs. fTR M to 



study equilibration phenomena in intermediate and high energy nucleus-nucleus collisions 
in a finite box with periodic boundary conditions. The semiclassical off-shell transport 
approach describes the virtual propagation of particles in the invariant mass squared M 2 
besides the conventional propagation in the mean-field potential given by the real part of 
the self energy. The imaginary part of the retarded self energy ImYT^p = — l/2Txp = 
—Pq(T co u(XP) + Tdecay{XP)) - apart from decay contributions to the width - is deter- 
mined by the collision integrals themselves and 'evaluated' within the transport approach 
dynamically. 

Our explicit calculations demonstrate that the off-shell propagation of nucleons has 
practically no sizeable effect on equilibration times especially at lower bombarding en- 
ergies; more importantly, the off-shell dynamics even lead to a slight increase of the 
equilibration time for kinetic equilibrium as noticed early by Danielewicz [|J (cf. also 
Ref. [^J). Furthermore, we have demonstrated that the off-shell HSD approach repro- 
duces the 'proper' spectral functions with respect to the statistical model (in case of a 
grand canonical ensemble) for nucleons and A's that in equilibrium are given analytically 
once the collisional width T co u(P, M) is known as a function of the 3-momentum P and 
invariant mass M. 



The authors like to acknowledge stimulating discussions with E.L. Bratkovskaya, C. 
Greiner and S. Leupold throughout this study. 
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Figure 1: The quadrupole moment in momentum space (|38|) for an infinite nuclear matter 
problem in a finite box with periodic boundary conditions characterized by bombarding 
energies of 0.1 A GeV and 1 A GeV at density p = p . The solid lines present the results 
for the off-shell calculations while the dashed lines are obtained in the on-shell limit. 
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Figure 2: The transverse mass distribution for an infinite nuclear matter problem at 
bombarding energies of 0.1 A GeV (upper part) and 1 A GeV (lower part) at density 
p = po- The solid lines present the results for the off-shell calculations while the dashed 
lines are obtained in the on-shell limit. 
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Figure 3: The differential distribution in mass for nucleons and A resonances at equilib- 
rium for an infinite nuclear matter problem at an initial bombarding energy of 1 A GeV 
at density p = po- The solid histograms present the results from the transport (box) cal- 
culations for the on-shell limit (upper part) and off-shell limit (lower part), respectively. 
The dashed lines ('thermo') are obtained from the statistical model at temperature T = 
97 MeV employing the spectral functions from the transport approach (see text). 
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Figure 4: Comparison of the differential distribution in mass for A resonances at equi- 
librium for the same infinite nuclear matter problem as in Fig. 3. The solid histogram 
presents the results from the transport calculations for the off-shell limit while the dashed 
histogram stems from the on-shell calculation. 
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